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Abstract 

Two-loop corrections with scalar and vector form factors are calculated for 
nuclear matter in the Walecka model. The on-shell form factors are derived 
from vertex corrections within the framework of the model and are highly 
damped at large spacelike momenta. The two-loop corrections are evaluated 
first by using the one-loop parameters and mean fields and then by refitting 
the total energy/baryon to empirical nuclear matter saturation properties. 
The modified two-loop corrections are significantly smaller than those com- 
puted with bare vertices. Contributions from the anomalous isoscalar form 
factor of the nucleon are included for the first time. The effects of the im- 
plicit density dependence of the form factors, which arise from the shift in 
the baryon mass, are also considered. Finally, necessary extensions of these 
calculations are discussed. 

( Figures are available by E-mail to HTANG@MPS.OHIO-STATE.EDU ) 
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I. INTRODUCTION 



The traditional theory of nuclear structure is based on the Schrodinger equation with a 
nucleon-nucleon (NN) potential that has its origin in meson exchange. In this nonrelativistic 
approach, one fits the NN potential to the empirical properties of the deuteron and to low- 
energy NN scattering data, and one then attempts to predict the behavior of many-nucleon 
systems. A natural and appealing generalization of this approach is to incorporate special 
relativity by using a relativistic quantum field theory with explicit meson and baryon degrees 
of freedom. These degrees of freedom are chosen because they are the most efficient for 
describing low- and medium-energy nuclear experiments. This generalization allows us to 
study interacting, relativistic, nuclear many-body systems, which future experiments will 
examine. The two basic questions are: what kind of field theories should we use, and how 
well can we describe nuclear systems using field theoretical models with hadrons, which are 
actually particles with internal structure. These are broad and hard questions that can be 
answered only after intensive and systematic investigation. 

Renormalizable relativistic quantum field theories with hadronic degrees of freedom, of- 
ten called quantum hadrodynamics or QHD, have been studied for some time P],0,[§ . At 
the level of the mean-field theory (MFT) and one-loop approximation ( "relativistic Hartree 
approximation" or RHA), these models can reproduce nuclear matter saturation and can 
realistically describe many bulk and single-particle properties of finite nuclei [@,[|f|]. The 
dynamical assumption behind renormalizability in QHD is that the quantum vacuum and 
the internal structure of the hadrons can be described in terms of hadronic degrees of free- 
dom alone. This assumption must ultimately break down at very short distances, and its 
limitations can be tested by explicit calculations. 

Some approximate QHD calculations beyond the one-loop level PJ^,0,U indicate large 
vacuum corrections, and even the validity of the one-loop vacuum contribution has been 
questioned |J. However, it is unlikely that we can calculate consistently beyond the mean- 
field level by introducing ad hoc procedures, such as including only the nucleons in the 
Fermi sea and simply throwing away Dirac-sea contributions. Indeed, it is already known 
that vacuum contributions are indispensable for maintaining the conservation of momentum 
and the electromagnetic current at the level of the random-phase approximation (RPA) 
[ jlfflnj . These results imply that we must develop practical and reliable techniques that go 
beyond the MFT and that include vacuum dynamics. 

A straightforward two-loop approximation for nuclear matter in the Walecka model was 
examined in Ref. ||, where large vacuum corrections were found. The loop expansion does 
not appear to be convergent or asymptotic in any sense. This is not a surprise if we notice 
that the two-loop corrections are essentially perturbative in the large couplings. An alterna- 
tive expansion in terms of meson loops (after integrating out the baryon fields) encounters 
the well-known ghost problems Even after the ghost poles are removed by Red- 

mond's "surgical" procedure ||12|| , putting aside for the moment its physical significance, the 
vacuum contributions are still big HJT3I]. These calculations suggest that vertex corrections 



(and short-range correlations) should be included to compute vacuum loops reliably in QHD 
theories |14 |. 



Since hadrons have internal substructure, one could argue for the use of nonlocal cou- 
plings and the introduction of ad hoc vertex form factors. This procedure was adopted in 
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Refs. |15j and . We remark that the introduction of form factors implies nonrenormaliz- 
ability of the theory from the outset, and thus the renormalization procedure in Refs. [HJ and 
|16| needs justification. (In other words, various subtractions were made in the computation 
of the energy that arise naturally only in a renormalizable theory.) This ad hoc procedure 
also says nothing about off-shell or density-dependent effects, and poses difficulties when 
one tries to construct approximations that obey the desired conservation laws. Moreover, 
although QCD in principle gives a complete description of nucleon structure, some part of 
the internal properties of the nucleon, particularly at large distances, must be equivalent to 
that provided by virtual hadron loops. These are the physical effects we intend to study in 
this paper. 

As pointed out recently by Milana ||14|| , a theory with baryons and vector mesons contains 
proper vertex functions that are highly damped at large spacelike momentum transfers, due 
to the contributions from virtual bremsstrahlung summed to all orders. This damping arises 
from the long-range (infrared) structure of the vertex and so should be calculable within 
the QHD framework. Since the QHD theory is renormalizable, the vertex function can be 
expressed in terms of the couplings and masses of the theory, with no ad hoc parameters. 
Although it may be impossible to achieve a truly quantitative description of the hadron 
structure within a QHD model, due to the complicated nature of the vertex functions, our 
goal is to understand at least the qualitative features of the vertex functions implied in a 
hadronic theory and to see how these features affect calculations with vacuum loops. 

the on-shell vector form factor was studied in a model with 
The behavior at large spacelike momentum transfer Q was 
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In a recent calculation 
baryons and vector mesons, 
determined by the leading logarithmic infrared behavior, which arises from the sum of all 
diagrams with vector ladders and crossed ladders across the hard vertex. At small Q, the 
lowest-order vertex correction was used, and the complete vertex function was constructed 
by interpolating smoothly between the low- and high-momentum transfer regimes. (See 
Ref. 



for a similar analysis with a somewhat different strategy.) 



In this paper we extend the method of Ref. [|17]] to the scalar-baryon vertex and apply the 
form factors so obtained to two-loop calculations in the Walecka model, with the lagrangian 
density 



C = ^{ifi-g y f-M + g s <t>)il) 
+ ±(«9 M 0«^0-m 2 2 ) 



where = d^V v — d v V^ and 8C contains the counterterms. We observe that the off-shell 
vertex functions should be used in a fully satisfactory calculation with loops. A full off- 
shell calculation is quite complicated, however, as one needs to know the off-shell behavior 
of the vertices at all spacelike momenta, as well as the modification of the vertices in the 
presence of valence nucleons at finite density. Calculations exploring these off-shell vertex 
functions are in progress [BJ|. Here, as a first step, we use an on-shell approximation, in 
which the off-shell vertex functions are replaced by their on-shell forms at zero density. This 
procedure is analogous to that used in Refs. []T3]| and [|I^| , where parametrized, on-shell form 
factors were used at the vertices, except that we use form factors obtained from within our 
model. Note also that the form factors used in Refs. and [|TBJ were chosen to have 
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simple momentum dependence, which allowed the two-loop energy to be calculated directly 
from the results of Ref. |J; in contrast, we allow for an arbitrary momentum dependence 
in the form factors and implement a totally new renormalization procedure for the two-loop 
contributions. Moreover, the anomalous vector form factor is included here in the two- 
loop results for the first time. Finally, we also estimate the effects of the implicit density 
dependence due to the effective nucleon mass M* that appears in the integrals that define 
the form factors. 

To include the form factors systematically, we apply the method of Freedman and McLer- 
ran [^] to nuclear matter in the Walecka model. One of the approximations that can be 
developed from this formalism is a two-loop approximation with one dressed vertex JTS 



Using this approximation, we calculate the corrections to the one-loop energy: first by using 
the coupling parameters and mean fields determined at the one-loop level and then by re- 
fitting the total energy to the empirical nuclear matter saturation properties. We find that 
the behavior of the two-loop corrections is significantly improved compared to that obtained 
with bare vertices 0. (This agrees with results obtained in calculations with ad hoc vertex 



functions fl5|,IE[|.) The contribution from the vector anomalous vertex is small compared to 



those from the vector charge and scalar vertices, but not negligible. The effects of the M* 
dependence and of the uncertainty in our interpolations for the form factors are modest com- 
pared to the overall size of the two-loop corrections, but are nevertheless significant on the 
scale of the nuclear matter binding energy. These new effects, which are usually neglected 
in studies using ad hoc form factors, are therefore important for a detailed description of the 
saturation properties of nuclear matter, because of the sensitive cancellations in the energy 
that occur near equilibrium density. It will certainly be necessary to extend our calculations 
to include the integrations over the off-shell vertex functions before any definitive statements 
can be made about nuclear matter saturation in this model. 

Some remarks are in order about our calculation of vacuum fluctuations involving the 
scalar field. Here we treat the scalar field that represents the a meson as "elementary", 
although it should, as generally believed, be considered as simulating the exchange of two 
correlated s-wave pions fl2l|,22,|23|,24l. Our goal in this work is to develop tools for calculating 
vacuum corrections in a model with an NN interaction that has short-range repulsion and 
mid-range attraction, as is empirically observed. The Walecka model is a simple one that 
satisfies this constraint. A more complete discussion of the role of chiral symmetry and the 
dynamical generation of the mid-range NN attraction, together with the implications for 
vacuum loops, should be considered as a necessary refinement of the work presented here 



Pj24]|. It is possible that a more detailed treatment could modify our results significantly, 
but it is also possible that a simple scalar field could remain an adequate approximation 
to correlated two-pion exchange even for the calculation of vacuum loops. We leave these 
extensions as topics for future study. 

The rest of this paper is organized as follows: In Sec. |I| we discuss the renormalization 
of the energy corrections in the two-loop approximation with one dressed vertex. The finite 
expressions that are used in our numerical calculations are generated. In Sec. [TT1| we present 
our approximations for the vertex functions and the method for determining the on-shell 
form factors. The results for the energy of nuclear matter with the vertex functions fixed 
at their free-space forms are given in Sec. [IV], and we study the effects of medium-modified 
vertices in Sec. [V] Section VI is a summary, and some technical details are included in an 
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Appendix. 



II. TWO-LOOP APPROXIMATION WITH VERTEX CORRECTIONS 

To include vertex functions in loop calculations of the nuclear matter energy density, 
it is convenient to apply the method of Freedman and McLerran ||2U|| , in which the ther- 
modynamic potential is constructed as a function of the full, connected propagators and 
proper vertices that satisfy the Schwinger-Dyson equations. By truncating the expansion 
of the thermodynamic potential and the Schwinger-Dyson equations appropriately, we can 
obtain various well-known approximations, such as the MFT, the RHA, the straightforward 
two- loop approximation ||, the relativistic RPA, etc. Details of these formal procedures are 



reported elsewhere |[13|J25|] . Here we merely quote the results for the two- loop approximation 
with one dressed vertex. 

Although the method of Freedman and McLerran provides for systematic truncation, it 
says nothing about the best way to approximate the exact thermodynamic potential in the 
case of strong couplings. It is also possible to truncate the expansion to include dressed 
vertices at both ends of the two-loop diagram. Here we shall take a conservative approach 
and include only one dressed vertex; if this gives adequate suppression of the vacuum loops, 
two dressed vertices will give even more. The question of which of these truncations is a 
better starting point can be answered only by including the next term in the expansion, and 
we leave this as a topic for future study. 

At the one-loop level (or RHA), the nuclear matter energy density can be written as |2| 

+ (^/ Fd3fc ^ + A£(AT), (2.1) 

where = (k 2 + M* 2 ) 1/2 and M* = M - g B cj) , with O the average scalar field. The 
spin-isospin degeneracy 7 = 4 for nuclear matter and 7 = 2 for neutron matter, and the 
mean vector field has been eliminated using 

since it is a constant of the motion. The one-loop vacuum correction ( "zero-point energy" ) 
is 

A£(M*) = - ^{M*Hn(M*/M) + M 3 (M-M*) - f M 2 (M - M*) 2 

+ f M(M - M*) 3 - f| (M - M*) 4 } , (2.3) 



where the renormalization conditions of Chin [2f| have been used, and we have assumed 
that 7 = 4 for the Dirac sea. 

The two-loop approximation with one dressed vertex obtained using the method of Freed- 
man and McLerran is most transparently described by Feynman diagrams. The two-loop 
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corrections to the RHA energy density can be drawn as in Fig. [I], where S denotes the 
renormalized proper baryon self-energy shown in Fig. [|, and A s (A v ) stands for the dressed 
scalar (vector) vertex. The square brackets indicate that the enclosed sub diagrams are 
renormalized by the inclusion of the appropriate subtractions. (We follow the diagrammatic 



conventions of Ref. |20fl 



The propagators in Fig. p] are as follows: The baryon Hartree propagator is 

G*(k) = G;(k) + GUk), (2.4) 

where 

Cm = - - * + M M * , , (2.5) 
k — M* + xe 

GUk) = -=;& + M*)5(k° - E* k )6(k F - |k|) , (2.6) 

and the momentum has been shifted to eliminate V , which is permissible for the evaluation 
of closed loops f| . The free (noninteracting) scalar and vector propagators are 

A°(g) = 2 \~. , (2.7) 

DlM=[-9^ + q -^)D°(q), (2.8) 

D°(q) = 2 \~. ■ (2.9) 

By separating the renormalized subdiagrams into unrenormalized parts and counterterm 
contributions (CTC), we obtain the more familiar diagrams in Fig. |3], where the CTC are 
not shown, since they can be determined easily from Fig. [I]. Evidently, an exact evaluation 
of these loops requires the knowledge of the off-shell proper vertices A s and A^, including 
their explicit density dependence. These off-shell vertices are currently under investigation 
I9| , but in this work, as a first approximation, we will replace them by their on-shell forms 



at zero density: 

ig a A a (q) = igsF s {-q 2 ) , (2.10) 
-ig y A>;(q) = -ig y [F vl (-gV 

+ iF v2 (-q 2 )a^q„] , (2.11) 

where a^ u = z[7 M ,7"]/2, and q^ is defined as the incoming momentum transfer at the 
vertex. Note that in the calculation of the energy, the on-shell form factors are needed only 
at spacelike momentum transfers, so we can use the results of Ref. |l7j . We emphasize that a 
similar on-shell approximation (with ad hoc form factors) is commonly used in essentially all 
calculations of nuclear matter properties (see, for example, Refs. JTJj and 0]); here, however, 



we will determine the on-shell momentum dependence within the context of our model. This 
also allows us to discuss some of the implicit density-dependent effects (contained in M*), 
which we examine in Sec. V. 
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With these considerations, we can translate the Feynman diagrams in Figs. |], |2], and |3| 
to obtain the renormalized two-loop correction 



+ l f4^ tT Pct^Wl - E - VEV > ( 2 - 12 ) 
7 Or „=i 

where "tr" denotes a trace over spin and isospin indices, and E ct contains the countert- 
erms for the baryon self-energy, as given in the Appendix. In Eq. ( |2.12| ), we have dis- 
carded the longitudinal term in the vector propagator, since its contribution vanishes 
by baryon current conservation, as can be easily verified in our approximation. The vacuum- 
expectation-value (VEV) subtraction is obtained by replacing all of the Hartree propagators 
G*(k) with noninteracting Feynman propagators G®(k) = Qt — M + ie)~ l and by omitting 
the quartic polynomial in <p . Note that some of the a n counterterms contain both one- 
loop [O(ft)] and two-loop [0(/i 2 )] contributions. The one-loop contribution comes from the 
renormalized meson polarizations in Fig. |I[ See Eq. (2.79) of Ref. H for the details of this 
separation. 

Due to our use of on-shell vertex functions in the two-loop integrals, the renormalization 
procedure is nonstandard and requires some discussion. In principle, the only unambiguous 
way to renormalize with vertex insertions is to use the off-shell vertices inside the integrals 
and include all required counterterms; this is clearly exceedingly difficult and motivates our 
simplified calculation using the on-shell vertices. We now observe that most of the effort 
in renormalizing the two-loop integrals with bare vertices involves the baryon propagators 
(in the self-energy and polarization loops). The noninteracting meson propagators simply 
follow along in the analysis, and although they determine which integrals diverge, they are 
otherwise innocuous. This suggests the following renormalization procedure with dressed 
vertices: 

1. We first carry out the renormalization of the vertex functions to arrive at finite, 
on-shell functions of the momentum transfer. This procedure is discussed in the 
next section, where we specify the detailed form of our model form factors. 

2. These form factors are then associated with the corresponding meson propagators 
(which carry the same momenta). The form factors then serve only to modify the 
momentum dependence of the meson propagators, producing so-called "M0ller 
potentials" |27|. 

3. The counterterms are then defined exactly as in the two- loop case with bare ver- 
tices ||, using the M0ller potentials in place of the noninteracting boson propa- 
gators. 

While this renormalization procedure may not be unique, it has the following advantages: 
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1. Finite results are obtained with the same number of counterterms as in the orig- 
inal two-loop calculation. All of the counterterms are constants that are defined 
by vacuum amplitudes, and all counterterm subtractions are local. 

2. All subintegrations for the meson polarizations and baryon self-energies satisfy 
the standard renormalization conditions || for any choice of vertex functions. 

3. The resulting energy density satisfies the same renormalization conditions as in 
the original two-loop calculation, for any choice of vertex functions. In particu- 
lar, the original results are obtained automatically when the model vertices are 
replaced by bare vertices. Note that, even with the form factors included, nested 
and overlapping subtractions remain divergent, although the overall subtractions 
are finite, the baryon self-energy counterterms. 

We remark that in a theory that postulates nonlocal vertices in the lagrangian, the 
definition of the renormalization counterterms is essentially arbitrary. In particular, since 
the overall subtractions to the energy are finite, there is no justification for stopping at 
0(<f)g), and one can simply remove all vacuum fluctuations by fiat. (That is, one can include 
an infinite number of counterterm subtractions.) To our knowledge, the only justification 
for stopping at a quartic polynomial is that the underlying theory is renormalizable, so that 
terms of 0(0^) and higher can never be removed. 

Following Ref. ||, we decompose the two-loop contribution into exchange, Lamb-shift, 
and vacuum fluctuation energies. Thus 

£ ^ = £ EX + £ls "I - £vf ; (2.13) 

where 

£ EX = \gl J 70j tr [Gl(k)Gl{k + q)] A°(g)F s (-g 2 ) 

-l^J^^^[Gmi,GUk + q )A^q)]D\ q ) , (2.14) 

£ L s = 9 2 s J tr [G;(k)GUk + q)] A°(g)F s (-g 2 ) + i j tr [E ct G* (*)] 

-9lJ^^tT[G;(kh,GUk + q)AmD°(q) , (2.15) 

SvF = \ 9 » I^Wf^ [G;{k)G;{k + q)] A °^(-<7 2 ) + 1 J (^4 tr PctG^)] 

" \ g " l-prr* ?^ tr raW^C* + <?) A ^)] D°(q) - E «„c - vev . 

2 J (2tt) 4 (2tt) 4 n=1 

(2.16) 

Here the divergent integrals have been regularized by writing them in r dimensions, with 
the limit r — ► 4 taken after the divergences have been removed. Note that these energy 
densities are implicitly functions of the baryon density p B and the mass parameter M*, with 
the latter to be determined by minimization of the full energy density. 
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The exchange contribution £ EX is finite. It is straightforward to work out the traces 
and arrive at the following form that is feasible for numerical integration using Gaussian 
quadrature: 



16tt 4 Jo 



dk 



k F pg2 

do 



+1 



dx 



+ 2g 2 F vl (P 2 



P 2 - 2M 
' P 2 + m 2 



*2 



+ 



' , 9N P 2 + AM* 2 
9sF s (P 2 ) p2M 2 

P + 771g 

6£ 2 M*P v2 (P 2 )P 2 



P 2 + m 2 



(2.17) 



where we have defined P 2 = 2(P^P* — M* 2 — kqx), with the dummy variables k = |k| and 
q = |q|. (Here and henceforth we take the spin-isospin degeneracy 7 = 4.) Note that since 
the momenta flowing through the meson propagators are spacelike, one does not encounter 
any poles. 

The Lamb-shift contribution £ LS can be written as 



d 4 k 



(2vr) 4 

S F a( 



tr[G*(fc)E P (fc)] 



M* 2 ) + M*S r 



*2\ 



-M 
-M* 2 



d 3 fc M 



(2tt)3 P* 
Ps(M*;p : 



-^ F -|k| 



+ M*S PB ( 

where the renormalized Feynman part of the baryon self-energy is 
S F (A;) = S FA (-A; 2 )+^ FB (-A ; 2 ) 



(2.18) 



(2.19) 



d T q 



g 2 s G;(k-q)A°(q)F s (-q 2 ) 
- glK{q)Gl{k - q)l»D\q) 



— E 



Ct ) 



and the specification of the counterterm coefficients 

S ct = -Cn(^-M)+ 7s o + M c 



(2.20) 



(2.21) 



is discussed in the Appendix. The scalar density of baryons is denoted by p s . 

By following the procedures in Ref . , the vacuum fluctuation energy £ VF can be evalu- 
ated by expanding the integrands in Eq. ( |2.16| ) in powers of (M* — M) using the algebraic 
identity 



G* F (k) = J2(M* - My[G° F (k)Y +1 + (M* - M) n+1 [G F (k)} n+1 G* F (k) 



(2.22) 



i=0 



Here n can be any positive integer or zero, and G F (k) is obtained from G*(k) by replacing 
M* with M. The zeroth-order term in the expansion [£ VP (M)] is cancelled by the VEV, 
and the counterterms a n are chosen as usual to cancel exactly the coefficients of the first 
four powers of (M* — M), which minimizes many-body forces p6|fl. The final result can 
be written as 
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x {gll l i[Gl(k)} 3 A^(q)D°(q) - g 2 [G F (£;)] 3 A (g)F s (-g 2 )) } 



+ \{M~ M*f /tJ^ tJ^ ^{[Gl{k - q)?Gl{k - q) 

x (^ 7 ,[G° F (A;)] 2 G;(A;)A^g)D^g)-^[G°(A;)] 2 G;(A;)A (g)F s (-g 2 ))} , 

(2.23) 

where is obtained from S K (fc) by replacing M* with M, and the renormalized vacuum 

scalar vertex at zero momentum transfer, 

A°(k)=A° A (-k 2 )+lfA B (-k 2 ) 

d T <? f „ 2 r^fO / 7_ _.\12 A 0/_.\ 771 / 2\ 



{g 2 M(k-q)] 2 ^(q)F s (-q 2 
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2 A^g)[G°(A;-g)] 2 7^)} 



7s 
9b 



(2.24) 



is discussed in the Appendix. Note that all integrals in Eq. ( |2.23j ) are finite, even for point 
vertices. 

After working out the traces, one can perform a Wick rotation to Euclidean momenta 
and evaluate the angular integrals analytically, leading to (£ VF = £ VF /M 4 ) 

(1 — m*) 5 y 



1 



dy 



2tx 2 Jo (y + l) 5 (y + m* 



dy 



(1 - m*) 5 y 



{S FA (y) [y 3 - 5(2 + m*)y 2 + 5(1 + 2m*)y - m* 

+ £° B (y) [(5 + m*)y 2 - 10(1 + m*)y + 1 + 5m*] y] 
{A° A (y) [(4 + m*)y 2 - 2(2 + 3m*)y + m* 



ds / dt 
Jo 



2tt 2 Jo " y (y + l) 4 (y + m* 2 ) 

-A° B (y) [y 2 - 2(3 + 2m*)y + 1 + 4m*] y) 

. -^{[t 3 -3{m* + l)t + m*][u 1 {8,t)A 1 {s) 

12o7r Jo Jo (t + l) 6 (t + m* z ) LL J 

-v 1 (s,t)B(s)}+ [(m* + 3)t -3m* - 1] [u 2 (s,t)A 2 (s) - v 2 (s,t)B(s)] } . 

(2.25) 

Here we have scaled all dimensional variables with the nucleon mass, set m* = M*/M, and 
defined 

g 2 F s (s) Ag 2 v F vl (s) 



A 1 (s) 



s + m? s + m 2 



g s 2 F s ( S ) , 2g 2 F vl (s) 
A 2 {s) = — : + 



B(s) 



s + m 2 s + m 2 
3ff 2 F v2 (s) 
s + m 2 



(2.26) 
(2.27) 
(2.28) 
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and 



Ui(s,t) 
u 2 {s,t) 



vi(s,t) 



v 2 (s,t) 



8st 
~2P ~ 

Ast, 
1 r 



m*C(s,t) 5 + m* 



1 



m 



*\2 



z 

t + m* 2 



s + t + l-Z) 



+ 



' 2(1 
s-t+3- 



- m 
1(1 



*A2 



C(s,*) 



*\2 



(s + t + l-Z) , 



Ast t — s + m 



+ 2 



*\2 



C(s,t) 



1 r 

Z 



t-s+3- Z 



[s + t + l-Z) 



4s* , 1X m*(t + s + m* 2 )^ / , 

Z?" + S + l) - 2fl- m -)' C(M) 



l 

2Z 



(5 + m*)(a + t + m* 2 - Z) + 9 + 2m* - 3m* 2 (s + t + l-Z) 



(2.29) 



C(s, t) = y/(s + t + m* 2 ) 2 - 4st - Z + — {s + t + 1) , 



(2.30) 



(2.31) 



(2.32) 
(2.33) 



with Z = y (s + 1 + l) 2 — 4st. We remark that the first two integrals in Eq. (|2.25 ) cor- 
respond precisely to the final two integrals in Eq. (2.93) of Ref. 0. The final integral in 
Eq. ( |2.25| ), however, was obtained here with a new renormalization procedure, and thus it 
cannot be directly compared to the quadrature in Ref. 0. 

The integrals in Eq. fl2.25| ) can be evaluated numerically using Gaussian quadrature. 
We found that splitting the integrations into two regions in each integral produced results 
that were accurate to better than 0.1% with a moderate number of points (~ 32 in each 
region). The largest uncertainty comes from the final integral in Eq. ( 2.25 ). We checked 
our computations by using two separate computer codes; all results for the nuclear matter 
energy presented below agreed to at least three digits. 



III. SCALAR AND VECTOR VERTEX FUNCTIONS 

In the preceding section, we renormalized the two-loop energy with the approximation 
that the fully off-shell vertices can be replaced by their on-shell forms at zero density, namely, 
Eqs. ( 2.101 ) and ( p. 11 ). Here we evaluate these form factors within the framework of our 



model. Since an exact calculation of the vertex function is impossible at present, we must 
make approximations. Our strategy is to include in the on-shell vertex functions the domi- 
nant physics that is accessible in a hadronic theory. 

We begin with the well-known fact that the proper vertex function in QED falls rapidly 
when the momentum entering on the photon line becomes large pH| , |2"9"|] . In particular, 
the asymptotic form for the on-shell vertex at large spacelike momenta q 2 < is 



11 



u{p b )A^u(p a 



u(p b )j>*u(p a ) 



x exp 



16tt 2 



In 2 



77T 



(3.1) 



where p\ 



Pb 



and m is an infrared regulator mass. In a theory with a massive 



neutral vector boson, the regulator mass m is replaced by the boson mass, and the electron 
charge e becomes the vector coupling g v . The physical origin of the strong damping is the 
large likelihood for virtual bremsstrahlung of soft vector bosons. In diagrammatic terms, the 
exponential arises from summing all ladders and crossed ladders involving the exchange of 
soft bosons across the single hard vertex, as shown in Fig. |j. We emphasize that although 
the momentum transfer to the vertex is large, the damping arises from the infrared structure 
of the theory, as the required factors of ln 2 (— q 2 ) are generated by loop momenta that are on 
the order of the vector meson mass. Thus it is reasonable to include this long-range vertex 
structure in a renormalizable theory containing hadron loops. Later work supports the 
assumption that non-leading logarithms appear only as multiplicative factors p0|j31| , |32|j33| . 

It is easy to show that the exponential damping in Eq. (|3.1|) is also reproduced by 
diagrams in which vector ladders and crossed ladders dress a single (hard) scalar vertex, 
at least to leading logarithmic order. Since the inclusion of higher-order vertex diagrams 
involving scalar meson exchanges produces only ultraviolet ( "hard" ) logarithms, as suggested 
by the work of Appelquist and Primack []3~4[| , these diagrams will not ruin the exponential 
damping |33|| . 

The general forms for the on-shell scalar and vector vertex functions at zero density are 
given in Eqs. ( p.lOj ) and (|2.11| ). Following Ref. JL7], we assume that the behavior of both 



vertices at small \q | is determined by the lowest-order vector vertex correction, that is, the 
middle diagram on the right-hand side of Fig. ^§ Then, by knowing the large \q 2 \ behavior 
of the form factors, we can interpolate smoothly between the small and large momentum- 
transfer regimes. Note that we need only spacelike (or Euclidean) momenta to compute the 
energy, so the form factors should be smooth, which makes the interpolation practical. 

The sum of the diagrams for either the scalar or the vector vertex in Fig. ^ results in 
a rapid suppression at large spacelike momentum transfers. As shown by Fishbane and 
Sullivan E9], the leading logarithmic asymptotic behavior for the on-shell vector charge 



form factor at zero density can be written as 



exp 



(3.2) 



where q = Pb — Pa and the superscript "(1)" indicates the lowest-order correction. Since the 
anomalous form factor is suppressed asymptotically by an additional factor of l/|g 2 |, we can 
write 



2MF v2 (-q' 



1 There is also a lowest-order vertex correction involving the exchange of a scalar meson. This 
additional diagram has a range similar to the diagram that we evaluate, so its inclusion would not 
significantly change the values that we obtain for the rms baryon radii. 
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exp 



(3.3) 



where the scale is set by the 0(g 2 ) anomalous moment, which is a conservative choice. 

At each order in g 2 , the dressed scalar vertex in Fig. |4] has the same denominator and 
the same leading power of q 2 in the numerator of the integrand as the dressed vector vertex. 
Since these features determine the asymptotic behavior p^ |, we conclude that the scalar 
and vector vertices behave similarly, and we can therefore write 



-Q 



exp 



-Q 



(3.4) 



which again holds for the leading logarithmic behavior. 

The lowest-order corrections to the on-shell vector charge and anomalous form factors 
are evaluated in Ref. |17|, and we simply quote the results here: 
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u 



1/2 



(3.7) 



Q 2 = —q 2 /M 2 , Q = y/Q 2 , and m v is written in units of the baryon mass M. 

The lowest-order correction to the on-shell scalar vertex can be calculated analogously. 
Thus 



F b {1) (-q 2 ) = Wlj-^k rG° F ( Pb - k)GU Pa - k)YD%(k) + CTC , 



(2k) 



(3.8) 



where the right-hand side is understood to be evaluated on shell {j> a = j> h = M), and the 
counterterm contribution (CTC) is determined by imposing the on-shell renormalization 
condition F^(0) = 0, which insures that the scalar-baryon coupling remains g s when the 
momentum transfer is zero. Straightforward manipulations similar to those for the vector 
vertex [|T7]] yield 



9t- 



16k 2 Jo 



du 



12S(u) 
Q 



In 



S(u) +uQ/2 
S{u) - uQ/2 



12u - 



Au(l — u + u 2 



u 2 + m 2 (l — u) 



+ 



2{l-u + u 2 ) + Q 2 {l-u + \u 2 ' 



QS(u) 



In 
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It is not hard to show |17| from Eqs. (|3.5|) and ( |3.9|) that for large Q 2 , 



Fil\Q 2 ), F^(Q 2 ) — -^^ 2 (Q 2 K) , (3-10) 



to leading order in logarithms. When combined with Eqs. (|3.2|) and (|3.4j ), this agrees with 
Eq. (O- 

As mentioned previously, we want to join the low- and high-momentum transfer regimes 
using a smooth interpolation. Evidently, the smaller the region to be interpolated, the more 
constrained the interpolation. On the other hand, the larger the high-momentum matching 
point Qo, the better the leading-logarithm asymptotic behavior. It is impractical, however, 
to take Qo so large that the non-leading logarithms are negligible compared to the leading 
logarithms. Thus, for a trade-off, we will approximate the non-leading logarithms by the 
large Q 2 behavior of FffiQ 2 ) and F S ^(Q 2 ). That is, we choose Qq large enough so that the 
values of Fj{ ] (Q 2 ) and F^(Q 2 ) can be accurately fitted by the function 

W(Q 2 ) = --f^ [\n 2 (Q 2 /ml) 

+ n HQ 2 /m 2 v ) + r 2 ] , (3.11) 

which means that all polynomials in 1/Q 2 are negligible. For Q > Q , we take the vertex 
form factors to be given by Eqs. ( |3.2| ), ( |3.3| ) and ( |3.4| ), with the lowest-order corrections 
replaced by Eq. ( |3.11| ). 



To minimize the uncertainty in the interpolations, one wants to choose Qo as small as 
is feasible, and we chose Qo = 5, as in Ref. [17]]. The parameters in Eq. (B . 1 1|) were 
determined by making a least-squares fit to the expressions in fl3.5| ) and (|3.9|)!i To test 
the sensitivity to the non-leading logarithmic behavior, we obtained one set of parameters 
for 5 < Q < 40, and then obtained another ("exact") set by letting Q 2 — > oo. (The 
exact T\ and r 2 can be determined analytically in principle, but the numerical evaluation 
is sufficiently accurate for our purposes.) Our first least-squares fit yields r\ rs —0.235 and 
r 2 ~ —1.60 for the scalar vertex, and r% ~ —2.88 and r 2 ~ 5.76 for the vector vertex, while 
the exact results are r\ = and r 2 = —3.35 for the scalar vertex, and r\ = —3.00 and 
r 2 = 6.66 for the vector vertex. (These values assume m v = 783 MeV and M = 939 MeV.) 
The resulting interpolations for the two choices of (as well as for ri = r 2 = 0) are 
nearly indistinguishable, which implies a very small sensitivity to the non-leading logarithmic 
behavior; for aesthetic reasons, we will use the "exact" parameters in all interpolations 
henceforth. 

In Ref. ||17|| , the last two terms in Eq. ( B.ll| ) were approximated by introducing a single 



parameter a such that the asymptotic behavior is exp[— (g^/167r 2 ) \r?(Q 2 /am 2 )}. For the 
scalar vertex, this is not as conservative as our present method. (That is, the method of 



2 We note that the onset of the asymptotic regime (Qo) an d the values of the parameters are 
determined by the behavior of the integrals in Eqs. ( |3.5| ) and ([O]). Thus these variables are all 
independent of the strength of the coupling g v . 



14 



Ref. leads to a scalar vertex function that decays more rapidly with increasing Q.) We 



therefore use Eqs. Q3.2Q , ( |3.3|) , and (|3.4j ) for Q > Q = 5, with the lowest-order corrections 
replaced by Eq. ( |3.11| ). We take the values of the form factors and their derivatives at Q = 0, 
as given by their lowest-order corrections, together with the values and the derivatives at Qq, 
determined from the asymptotic forms, as four input parameters to specify the interpolation 
functions. 

In Ref. [F7|, the interpolating function was taken as 

fu{Q 2 ) = aexp(-6Q 2 ) + c/(l + Q 2 ) + d , (3.12) 

where a, b, c, and d are to be determined. This interpolating function contains both gaussian 
and monopole terms, so we will call it a "mixed" interpolation for convenience of description. 
As noticed in Ref. |T7|], for some values of g 2 , two different solutions for the parameters a, b, 
c, and d can be found, while for other values, unique solutions exist. Whereas this behavior 
can furnish an estimate of the uncertainty in the interpolation, it also results in a nuclear 
matter energy that is a discontinuous function of g 2 . This is inconvenient when one attempts 
to refit the couplings to reproduce nuclear matter saturation. 

Thus we use instead the following polynomial interpolating function, 

/p(Q 2 )=a/(l + Q 2 ) 3 + 6/(l + Q 2 ) 2 

+ c/(l + Q 2 ) +d , (3.13) 

which yields a unique solution for a, b, c, and d at any g 2 . To estimate the uncertainties in 
the interpolation and their effects on the energy density of nuclear matter, we also use the 
gaussian interpolating function 

f G (Q 2 ) = aexp(-Q 2 ) + &exp(-Q 2 /5) 

+ cexp(-Q 2 /10) + d , (3.14) 

which also yields a unique solution for a, b, c, and d at any g 2 . Since F v2 (g 2 ) is suppressed 
asymptotically by a factor of l/\q 2 \ relative to F vl (g 2 ), we take its corresponding interpolat- 
ing functions to be those of Eqs. ([TLj) and (gTLj ) multiplied by 2MFi 2 1) (0)/(l + Q 2 ^ 



In Fig. H] we show the scalar form factor and the vector charge form factor using the 
mixed, polynomial, and gaussian interpolating functions for RHA couplings. The vector 
anomalous form factors obtained from Eqs. ( 3.1 3| ) and ( |3.14 ) with a further multiplicative 



factor of 2MF V ( 2 1) (0)/(1 + Q 2 ) are shown in Fig. |6|(a). The uncertainty envelopes generated 
by the polynomial and gaussian interpolations are combined in Fig. |6|(b). In Fig. 0(a) we 
show the vector charge form factor and the corresponding monopole form factor [F m (Q 2 ) = 
A 2 /(Q 2 + A 2 )] that decays similarly, while in Fig. 0(b), we show similar results for a dipole 
form factor [F&{Q 2 ) = A 4 /(Q 2 + A 2 ) 2 ]. The parameters a, b, c, and d for the preceding 
interpolations are listed in Table [| All values assume g 2 = 102.8, m v = 783 MeV, and 
M = 939 MeV. 

The preceding figures indicate that the decay of the form factors as a function of Q 
depends significantly on the interpolation functions. Nevertheless, we have verified that 
the polynomial and gaussian forms provide a reasonable envelope on the uncertainty intro- 
duced by the interpolation. We will therefore use the polynomial and gaussian functions 
to investigate the sensitivity of the nuclear matter binding energy to the decay of the form 
factors. 
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IV. RESULTS 



We now present numerical results for the vertex-corrected two-loop energy computed with 
the equations in Sec. H an d vertices in Sec. |T|. We start with the results of a perturbative 



calculation, in which the RHA couplings and self-consistent mass M* are used to compute 
the two- loop corrections. For convenience of comparison with the results of Ref. |J, we use 
the couplings gl = 54.3 and g% = 102.8, and masses m s = 458 MeV and m v = 783 MeV, 
which produce nuclear matter equilibrium at k F = 1.30 fm -1 with a binding energy/nucleon 
of 15.75 MeV. Later we present results obtained by refitting the parameters so that the total 
energy, including the two- loop contributions and the new self-consistent M*, reproduces the 
preceding equilibrium properties of nuclear matter. Note that all vertex functions used in 
this section are fixed at their free-space parametrizations, which follow from the formulas 



in Sec. |T|. Calculations both with and without the anomalous form factor F v2 will be 
considered. 

The results for the perturbative two-loop energy/baryon of nuclear matter as a function 
of density, with both the polynomial and gaussian interpolations for the vertices, are shown 
in Fig. |8]. As a check of our renormalization procedure and numerical calculations, we 
reproduced the two-loop results with bare vertices as shown in Ref. || by setting F s = 
F vl = 1 and F v2 = O.i By comparing the results of Fig. |S] with those in Fig. 4 of Ref. ||, we 
observe that the form factors suppress the vacuum contributions considerably. Moreover, 
the difference between the two interpolations gives an estimate of the uncertainty in the 
model vertex functions, which produces large changes on the scale of the nuclear binding 
energy. Thus the energy density of nuclear matter in this model is sensitive to the behavior 
of the form factors at intermediate momenta, at least in the two-loop approximation. 

The separate contributions from the scalar, the vector charge (F Y i), and the vector 
anomalous {F v2 ) parts are shown in Fig. ^. Evidently, there are strong cancellations be- 
tween the vector charge and scalar pieces. Notice also that the anomalous contributions are 
generally only 5%-15% as large as the scalar and vector charge contributions. Neverthe- 
less, the anomalous terms are comparable to the binding energy of nuclear matter and are 
relevant in a calculation of nuclear matter saturation because of the sensitive cancellations 
between the vector charge and scalar contributions. 

The separate contributions to the two-loop correction from the vacuum fluctuation, the 



exchange, and the Lamb-shift terms are shown in Fig. [L0] with F v2 omitted and in Fig. [LI 
with F v2 included. One observes that the Lamb-shift energy suffers less suppression from the 
form factors than the vacuum fluctuation energy does. When F v2 is included, the vacuum 
fluctuation and Lamb-shift energies become more negative. 

The two-loop energy with vertex modifications £^ 1 ' can be minimized with respect 

to M* at each density, and the parameters can be adjusted to reproduce the equilibrium 
properties of nuclear matter. The vector meson mass is held fixed at its empirical value of 
783 MeV (for simplicity), but the scalar meson mass must be increased, with the size of the 



3 We can also reproduce the results in Ref. [l5| by setting F v2 = and choosing F s and F v i equal 
to the square of the monopole form factors therein. 
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increase determined by how big the two-loop perturbative result deviates from the RHA 
result. (If too small a value is chosen for the scalar mass, the resulting interaction is too 
attractive, and nuclear matter always saturates with too large a binding energy.) 

We have found a number of parameter sets that reproduce the equilibrium properties 
of nuclear matter. A representative sample is listed in Table [TT], where the compressibility 
K~ l is also shown.! Fig. O shows the corresponding refitted results with and without the 



anomalous form factor and with the polynomial or the gaussian interpolation. In general, 
the compressibility remains large in this approximation and increases slightly from the value 
in the RHA; the form factors with the polynomial interpolation produce more suppression 
of the vacuum contributions (compared to the gaussian interpolation), but give a higher 
compressibility. The corresponding self-consistent nucleon mass is shown as a function of 
density in Fig. [T^. Here M* is obtained by minimizing the full two-loop energy density 
£W(M*) + £ ( 2 )(M*) with respect to this parameter. The two-loop contributions generally 
reduce the value of M* . 

The RHA and two-loop contributions to the total refitted energy with F v2 excluded are 
shown in Fig. |TJ], where we also show the vacuum fluctuation, the exchange, and the Lamb- 
shift energies. The RHA and two-loop contributions to the total refitted energy with F v2 
included are shown in Fig. where the scalar, vector charge, and vector anomalous parts 
in the two-loop correction are also shown. 

To refit the energy/baryon to nuclear- matter equilibrium, the scalar mass must be in- 
creased significantly, which suggests that the results in this approximation are sensitive to 
the precise shape of the form factors. Nevertheless, the two-loop calculation with vertex 
modifications is a better approximation than the conventional two-loop calculation with 
bare vertices, since the two-loop corrections are smaller. For quantitative comparisons to 
nuclear matter properties, one must know the form factors more accurately, especially in 
the intermediate momentum transfer region. Thus our calculations imply that it is inap- 
propriate to draw conclusions from calculations using some particular ad hoc choice of form 



factors, as in Refs. |15j and [116]. 



V. RESULTS WITH MEDIUM-MODIFIED VERTICES 

One advantage of having an explicit model for the meson-baryon vertex functions is that 
we can investigate how these functions change in the nuclear medium. For example, the 
baryon mass changes from M to M*, and by making this change in the baryon propagators 
that appear in the vertex loops, the vertex functions acquire an implicit density dependence. 
This is the effect we will consider in this section. We emphasize that this is a first step in the 
study of density-dependent vertex modifications in this model, since we are neglecting the 
valence-nucleon contributions to the vertex loops, as well as additional form factor functions 
that can arise at finite density. Moreover, as with the calculations in the preceding sections, 
a quantitative study requires the full off-shell vertex, as well as the inclusion of pions to more 



4 As is well known, the compressibility can be reduced, if desired, by including nonlinear scalar 
meson self- interactions. 
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accurately describe the long-range vertex structure. We leave these additional modifications 
as important topics for future investigations. Although the present calculation is just a 
first step, it illustrates numerous issues that must be dealt with in any microscopic model 
of the meson-baryon vertices; many of these issues are simply omitted by assumption in 
conventional calculations based on ad hoc form factors. 

By replacing M with M* in the vertex functions, the mass of the virtual intermediate 
state is reduced, and thus the radius of the dressed nucleon increases with increasing density. 
This is the primary modification we study here. We emphasize that the leading logarithmic 
behavior at large spacelike Q 2 is unchanged. However, since we determine the non-leading 
behavior from an examination of the lowest-order loop diagrams, this behavior also changes 
at finite density; these effects are very small and are incorporated here just for consistency. 
We will follow the general strategy described earlier: include effects that can be calculated 
reliably in our simple model, and treat effects that cannot yet be calculated as conservatively 
as possible. 

It is a straightforward exercise to generalize the results in Sec. || to incorporate the 
modified baryon mass, and we simply quote the results here: 
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M*/M, and m v is written in units of the baryon mass M. 



We note the important point that the prefactor 2M is used solely to make the anomalous 
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contribution dimensionless (any mass could be used), and thus it does not change to M* at 
finite density. 

Here we have renormalized the form factors using the same renormalization conditions 
as in Sec. [HI[ A straightforward numerical evaluation of the preceding integrals shows 
that they do not vanish at Q 2 = 0, and thus the vector charge and scalar couplings are 
renormalized at finite density. This is especially surprising for the vector term, since the 
conservation of the baryon current implies that there is no charge renormalization. The 
resolution of this dilemma is simple. In a calculation that includes the full off-shell vertex 
inside the two-loop diagrams, one should also include the appropriate self-energy insertions 
on the baryon lines to maintain the required Ward identities. Thus, when we "extract" the 
vertex and renormalize it on mass-shell, using the procedures discussed in Sec. U, we should 
also extract the self-energies, which now appear on external baryon lines. In free space, 
these self-energy insertions vanish by construction, due to our choice of mass and wave- 
function counterterms, but when M is replaced by M*, one finds a finite wave- function 
renormalization. It is straightforward to show |35] that the contribution from this wave- 
function renormalization precisely cancels the shift from F*^(0), and the baryon number 
remains unchanged. Thus we will set the vector charge form factor F* x (0) = 1 in performing 
the interpolations discussed below. 

The self-energy insertions also produce a finite correction to the baryon mass, which is 
exactly what is needed to modify the mass in the external spinor to include the correction 
from the exchange self-energy. Since we will only consider perturbative calculations about 
the RHA results, in which we neglect all exchange corrections to the baryon mass (by 
definition), we will neglect this mass shift as well. 

In contrast, the scalar vertex is not protected by the Ward identity. The finite wave- 
function renormalization does not cancel exactly against F*W(0), and there is a finite shift in 
the strength of the vertex. However, while we believe that the lowest-order vertex correction 
gives a reasonable estimate of the form-factor radius (since it incorporates the intermediate 
state with the lowest mass), there is no reason to expect that it is reliable for the renormal- 
ization of the scalar coupling strength. This renormalization is a completely dynamical effect 
that is likely to depend significantly on valence-nucleon contributions, scalar meson exchange 
diagrams, vacuum polarization, etc., all of which are neglected here. We will therefore take 
a conservative approach and leave the magnitude of the scalar strength unrenormalized at 
finite density; we postpone the difficult problem of obtaining a reasonable estimate for the 
finite renormalization to a future investigation. 

It is also straightforward to obtain expressions for the medium-modified rms radii of 
the nucleon, but we will not present the formulas here. (These are most easily obtained 
by differentiating the Feynman parameter integrals with respect to Q 2 , rather than by 
differentiating the expressions given above.) One finds that the vector charge and scalar 
mean-square radii scale approximately as 1/M*. Because of the overall factor of m* in 
2MF* 2 , however, the anomalous moment (and its radius) are insensitive to the value of M*, 
and this form factor is essentially the same as at zero density. 

The parametrizations of the vertex functions are performed as described in Sec. PTT| , with 
a different parameter set at each value of M*. One finds that the values of r\ [see Eq. ( |3.11D 1 
are unchanged at finite density, while the values of r2 suffer a (small) change that is exactly 
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the same for the scalar and vector vertex. Fig. [16] shows the vector charge form factor 
at M* = 0.7 M, where we have included the finite wave-function renormalization, so that 

f; 1 (o) = i. 

If one faithfully carries out the renormalization procedure defined in Sec. [II] using the 
medium-modified vertices, one discovers that the resulting expressions for the two-loop en- 
ergy are no longer finite. This occurs because the two-loop integrals contain the modified 
vertices, while the subtraction terms are defined by vacuum amplitudes and use unmodified 
vertices; thus, certain divergences no longer cancel. The problems can be traced to the nested 
divergences, which require the renormalization of the self-energy and vertex functions, as 
discussed in the Appendix. (The overlapping subtractions remain unchanged, and the over- 
all subtractions are always chosen to remove the first four powers of <f> .) In principle, the 
only way to determine the required new subtractions is to use the full off-shell vertices inside 
the loops. Within the context of the on-shell approximation for the form factors, the best 
that can be done is to use medium-modified vertices for both the two-loop integrals and the 
corresponding subtractions. This is the procedure we follow here. The result is that the en- 
ergy density is calculated by inserting the density- dependent vertex functions, parametrized 
as discussed earlier, into the expressions for the renormalized two-loop contributions given 
in Sec. [TJ. 

In Fig. |T7|, we show the results of a perturbative calculation (RHA parameters and M*) 
using the medium-dependent vertices with both the polynomial and gaussian interpolations. 
A comparison with Fig. |8| shows that the two-loop corrections are smaller here, as expected 
from the larger radii (and consequent more rapid decay in momentum space) of the F* t 
and F* vertices. The differences between the results with the two different interpolations 
is again an estimate of the uncertainty in the form factors at intermediate momenta; these 
differences remain significant on the scale of the nuclear matter binding energy. Moreover, 
it is clear that including the density dependence of the vertices also produces effects that are 
significant on this scale. This suggests that for quantitative calculations of nuclear matter 
saturation, it is necessary to include this density dependence, and one can question the 
accuracy of results using vertices that are simply assumed to be density independent. 



VI. SUMMARY 

We have calculated the two-loop energy with one dressed vertex for nuclear matter in 
the Walecka model, using form factors determined in the framework of the model. With an 
on-shell approximation for the vertex, we determine the leading behavior at large spacelike 
momentum transfer from the sum of an infinite set of ladder and crossed- ladder diagrams, 
and we use the lowest-order vertex correction to specify the behavior at small momentum 
transfer. Two different interpolation functions are used to join the two momentum-transfer 
regimes. Whereas the resulting vertex functions are not determined uniquely at intermediate 
momentum transfers and are known only to within an error band, they are specified by the 
underlying couplings and masses in the model without the introduction of ad hoc parameters. 

We find that the two-loop corrections to the one-loop energy are considerably smaller 
than those computed earlier with bare vertices. While these corrections produce significant 
changes in the nuclear binding energy, one can still fit the empirical equilibrium properties of 
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nuclear matter with a reasonable adjustment of the model parameters. (The most dramatic 
adjustment is in the scalar meson mass, which must be increased by roughly several hundred 
MeV.) Since these results were obtained using only one dressed vertex in the loop, it is likely 
that calculations using two dressed vertices will find even smaller corrections. 

We also included the contributions from the anomalous vector vertex for the first time. 
Although its contribution is small compared to those of the vector charge vertex or the scalar 
vertex, it is not negligible on the scale of the nuclear matter binding energy. Moreover, the 
implicit density dependence in the vertices and the uncertainty in our interpolations also 
produce variations that are important on the scale of the binding energy. These new effects 
should be present in any realistic model of the meson-baryon vertex functions and must be 
included before a quantitative calculation of nuclear matter saturation can be performed. 

Although the present calculations are the first step in the inclusion of dynamically gen- 
erated vertex functions in relativistic nuclear matter calculations, and many improvements 
and refinements must be made, there are several features of our analysis that we believe to be 
under control. First, the leading logarithmic behavior of the vertices at large spacelike mo- 
mentum transfer is included correctly. Second, the low-momentum-transfer behavior should 
be reproduced reasonably well (for a model with only heavy mesons), since we include the 
intermediate state with the smallest mass. Third, our two different interpolation functions 
give a realistic "error envelope" for the uncertainty of the vertex functions at intermediate 
momentum transfer and for the resulting uncertainty in the nuclear matter binding energy. 
Fourth, the calculated isoscalar anomalous moment has a magnitude that is consistent with 
empirical values, and thus it provides a meaningful estimate of the size of the isoscalar 
anomalous contributions to nuclear matter saturation. Finally, the mean-square radii of the 
vector charge and scalar vertices increase with density roughly as the inverse of the baryon 
effective mass, which should give sensible results for the impact of this density dependence 
on the nuclear matter calculation. 

There are several improvements that must be made in these calculations before any defini- 
tive conclusions can be drawn. First, one must compute the off-shell vertices and include 
them inside the loop integrals; this will also require a computation of the corresponding 
self-energies, which must be included to maintain the conservation of the baryon current 
||19|| . Second, the contributions from the density-dependent (valence) parts of the baryon 
propagators must be included in the computation of the vertex functions; this will introduce 
(in principle) new vertex functions that arise at finite density. Third, pions must be in- 
cluded to describe the long-range vertex structure more accurately. Fourth, the description 
of the vertex functions at intermediate momentum transfer should be improved; this may 
be possible by examining the dispersion relations that determine these functions. Finally, 
one must investigate the truncation procedure used to define the present approximation to 
decide if it is more accurate to dress one or both of the vertices in the two-loop terms. All of 
these improvements can be studied systematically within the QHD framework and provide 
topics for future work on this problem. 

In a larger context, regardless of which degrees of freedom one believes are the most 
appropriate, it is necessary to have reliable models of the off-shell and density- dependent 
behavior of the meson-baryon vertices before accurate relativistic calculations of nuclear 
matter properties can be made. The construction of such models presents a formidable 
challenge to the practitioners of relativistic nuclear many-body theory. 
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APPENDIX: BARYON SELF-ENERGY AND BARYON SCALAR VERTEX 



Here we derive the renormalized Feynman parts of the baryon self-energy and baryon- 
scalar vertex that appear in Eqs. ( |2.20| ) and ( |2.24| ). The baryon self-energy is shown graphi- 
cally in Fig. H The baryon-scalar vertex is shown in Fig. [18], with the counterterm subtrac- 
tion omitted. Here we will be interested only in the case of zero momentum transfer. Our 
renormalization procedure is similar to that in Ref . [0] . 

It follows from Eq. (|2T2q) that 



S F (fc) — S FA ( k ) ~f" ^E FB ( k 
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The convergence properties of these integrals depend on the behavior of the form factors; 
for bare vertices, the integrals are logarithmically divergent. To renormalize them, we write 
them in terms of Euclidean momenta and then explicitly subtract the integrands using 
counterterms chosen to reproduce the appropriate renormalization conditions. (This method 
is discussed in Refs. and |jl3| .) One can verify that, for bare vertices, this procedure 
reproduces the results obtained with the more conventional dimensional regularization ||, 
and it allows us to extend the renormalization to include our model vertex functions. 

After making a Wick rotation to Euclidean space (po ip%) and performing the angular 
integrals, we can write Eqs. (|A"2|) and (|A3|) as 
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where A 1 (t),A 2 (t), and B(t) are defined in Eqs. (|2~2^| ), (|2^7| ), and (|2~2^| ), and 
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Here we have set s = k E and t = q E 2 , with the subscript "E" denoting Euclidean. 
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Similarly, from Eq. ( |2.24|) , we have 

A (k)=A° A (-k 2 )+k'A B (-k 2 ) , 
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51), and ( |A9| ) are defined by imposing the 
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After inserting the expresions for the self-energy and vertex [Eqs. ( |AT| ) and (|A8|) 1 into 
Eqs. ( |A13|) — (p\15|) , it is straightforward to find the counterterm contributions and to make 
the necessary subtractions. After some tedious algebra, the final results are 
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where we have scaled all dimensional variables with the nucleon mass M and defined r] t = 
yl +4/t and m* = M*/M. 

All of the preceding integrals can be evaluated using Gaussian quadrature, although 
care must be used due to cancellations at large t. It is easily verified that the integrals that 
contain F s and F vl are all finite even after setting F s — 1 and _F vl = 1, while the integrals that 
contain F v2 are also finite when one realizes that in any realistic calculation, F V 2(s) / F v i(s) 
is suppressed by at least a factor of 1/s at large s. We have verified numerically that our 
results with F s = F vl = 1 and F v2 = are the same as those of Ref. 0. 
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FIGURES 



FIG. 1. The two- loop corrections to the RHA energy density. The solid, dashed, and wavy lines 
represent the baryon Hartree propagator, the free scalar propagator, and the free vector propagator, 
respectively. The subtractions for the overall divergences and for the vacuum expectation value of 
the energy are not shown. 

FIG. 2. The proper baryon self-energy. The diagrammatic notation is the same as in Fig. |]. 

FIG. 3. The unrenormalized two-loop diagrams. 

FIG. 4. Diagrammatic expansion of the proper scalar and vector vertices. 

FIG. 5. The on-shell scalar form factor (a) and vector charge form factor (b) with the mixed, 
polynomial, and gaussian interpolations. 

FIG. 6. The on-shell vector anomalous form factor (a) and the uncertainty envelopes generated 
by the gaussian and polynomial interpolations (b). 

FIG. 7. Comparison of the on-shell vector charge form factor of Fig. ||(b) with some similar 
monopole and dipole form factors. 

FIG. 8. The perturbative two-loop energy with vertex corrections using the polynomial and 
gaussian interpolations. The dotted (dashed) curves are obtained with F v 2 excluded (included). 
The solid curve is the RHA energy. Both calculations use RHA parameters and m v = 783 MeV. 
The equilibrium density p corresponds to k F = 1.30 fm — . 

FIG. 9. The separate contributions to the total two-loop correction from the scalar, vector 
charge, and vector anomalous parts. The parameters are the same as those in Fig. ||. 

FIG. 10. The separate contributions to the total two-loop correction from the vacuum fluctu- 
ation, the exchange, and the Lamb-shift terms, excluding F v 2- The parameters are the same as 
those in Fig. ||. 

FIG. 11. The separate contributions to the total two- loop correction from the vacuum fluctu- 
ation, the exchange, and the Lamb-shift terms, including F v 2- The parameters are the same as 
those in Fig. 

FIG. 12. The refitted two-loop energy with vertex modifications obtained for the polynomial 
and gaussian interpolations, respectively. The parameter sets correspond to those in Table O. 
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FIG. 13. The effective nucleon mass for the refitted two-loop approximation with vertex mod- 
ifications obtained for the polynomial and gaussian interpolations, respectively. 

FIG. 14. The RHA and the two-loop contributions to the refitted energy density with F v 2 
excluded. The vacuum fluctuation, the exchange, and the Lamb-shift parts of the two-loop contri- 
bution are also shown. 

FIG. 15. The RHA and the two-loop contributions to the refitted energy density with F v 2 
included. The scalar, vector charge, and vector anomalous parts of the two-loop contribution are 
also shown. 

FIG. 16. The vector charge form factor F v i at M* = 0.7M. The solid (dashed) line is obtained 
with polynomial (gaussian) interpolation. The dotted and dot-dashed lines are the corresponding 
vacuum form factors, which are the same as those in Fig. |5](b). 

FIG. 17. The perturbative two-loop energy with density-dependent vertices. The couplings are 
the same as those in Fig. ||| The dotted (dashed) curves are obtained with F v 2 excluded (included). 
The solid curve is the RHA energy. 

FIG. 18. The baryon-scalar vertex that corresponds to Eq. ( p. 24 ). 
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TABLES 



TABLE I. Interpolation parameter sets. 





a 


b 


c 




d 




mixed F s 


0.747185 


0.457293 


2.60765 xlO" 


- 1 


-7.94998 xlO" 


Q 

-o 


F vl 


0.792431 


0.344247 


2.12528 xl(T 


-i 


-4.95969 xlO" 


-3 


F s 


-1.30867 


2.21488 


9.87004 xHT 


-2 


-4.91063 xlO" 


-3 


poly. F vl 


-1.45088 


2.38433 


6.92982 xlO" 


-2 


-2.75043 xlO' 


-3 


F v2 


-2.94795 


3.98240 


-3.35013 xlO" 


-2 


-9.41037 xlO" 


-4 


F s 


0.498343 


0.541225 


-4.14082 xlO" 


-2 


1.83978 xlO" 


-3 


gauss. F v i 


0.348549 


0.720836 


-7.39583 xlO' 


-2 


4.57337 xlO- 


-3 


F v2 


-1.43842 


2.83889 


-4.19250 xlO" 


-1 


1.87796 xlO" 


-2 



TABLE II. Two- loop refitted parameter sets. 



set 


interp. 


F v2 




9l 


m s (MeV) 




A 


poly. 


No 


63.6 


130.8 


458 


600 


B 


poly. 


Yes 


102.6 


109.6 


600 


524 


C 


gauss. 


No 


159.4 


79.4 


800 


437 


D 


gauss. 


Yes 


233.6 


106.2 


900 


494 
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